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We present a detailed analysis of the electronic and optical properties of two-electron quantum dots 
with a two-dimensional Gaussian confinement potential. We study the effects of Coulomb impurities 
and the possibility of manipulate the entanglement of the electrons by controlling the confinement 
potential parameters. The degree of entanglement becomes highly modulated by both the location 
and charge screening of the impurity atom, resulting two regimes: one of low entanglement and 
other of high entanglement, with both of them mainly determined by the magnitude of the charge. 
It is shown that the magnitude of the oscillator strength of the system could provide an indication of 
the presence and characteristics of impurities that could largely influence the degree of entanglement 
of the system. 
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I. INTRODUCTION 

Semiconductor quantum dots (QDs) are excellent can- 
didates for realizing qubits for quantum information pro- 
cessing because of the potential for excellent manipu- 
lability and scalability. In contrast to real atoms and 
molecules, in QDs the electronic and optical properties 
are highly tunable. Tremendous advances in semiconduc- 
tor technology allow the preparation of complex struc- 
tures and give the possibility to the experimentalists to 
have a great control on the parameters that define the 
electrical and optical properties of these systems^. 

It is known that the presence of impurity centers has 
a great influence on the optical and electronic properties 
of nanostructured materials. Since the pioneering work 
of Bastard^ many authors have investigated the effects of 
impurities on different properties of artificial atoms and 
molecules. A recent work^ studied the effects of hav- 
ing unintentional charged impurities in two-electron lat- 
erally coupled two-dimensional double quantum-dot sys- 
tems. They analyzed the effects of quenched random- 
charged impurities on the singlet-triplet exchange cou- 
pling in two-electron double quantum-dots. Although 
there is an enormous interest in applying these systems 
in quantum information technologies, there are few works 
trying to quantify the effect of charged impurities on this 
kind of tasks. The existence of unintentional impurities, 
which are always present in nanostructured devices, af- 
fects seriously the possibility of using these devices as 
quantum bits. Although the distribution and concentra- 
tion of impurities in these systems result unknown pa- 
rameters, there are some recent works that propose the 
possibility of experimentally control these issues^—. Im- 
purity doping in semiconductor materials is considered 
as a useful technology that has been exploited to control 
optical and electronic properties in different nanodevices. 

The entanglement, which is one of the most curious 
phenomena in quantum mechanics, is being considered in 
recent years as a physical resource that can be used for 



quantum information processes as teleportation of quan- 
tum states^—. There exists the possibility of manipu- 
late the amount of entanglement in a QD molecule by 
controlling the nanostructure parameters that define the 
nanodevice. Zunger and He 11 studied the effect of inter- 
dot distance and asymmetry on the spatial entanglement 
of two-electron coupled quantum dots. They showed 
that the asymmetry in these systems significally lowers 
the degree of entanglement of the two electrons. Two- 
electron entanglement of different quantum dots atoms 
and molecules have been studied by several authors in 
the last decade^ - —. The presence of nearby charged 
impurities, as Das Sarma and Nguyen show, have an im- 
portant effect on the singlet-triplet coupling, with un- 
wanted impact in quantum information tasks. One of 
the main goals of this paper is the calculation of the spa- 
tial entanglemen t 10 i 13 ' 15 of the two electrons in a double 
QD molecule in presence of charged impurities. Experi- 
mentally, results very difficult to mesasure the amount of 
entanglement of two electron in a coupled QD directly. 
There exist several techniques that allows one to mea- 
sure the possibility of double occupation^ and the opti- 
cal properties such as the dipole transition, the oscilla- 
tor strength and the photoionization cross sectionSIr— of 
these systems. If we know the relationship between these 
quantities and the degree of entanglement, we can have 
information about the amount of entanglement of our 
system. The possibility of using this information in order 
to design nanodevices according to the level of entangle- 
ment desired results quite difficult because the positions 
and the strength of the impurities are unknown. Despite 
this, there are some recent experiments which show the 
mechanism of dopant incorporation and how the incor- 
porations of impurity deflects can be controlled-£i22. 

The aim of this work is to present a detailed analysis of 
the electronic and optical properties of a two-dimensional 
two-electron coupled quantum dot and the effect of im- 
purties. In particular we show that the entanglement of 
the electrons is strongly modulated by the position and 
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charge of the impurity. We also show that optical mea- 
surements would allow to obtain information about the 
effect of the impurity in these kind of devices. The pa- 
per is organized as follows. In Section II, we introduce 
the model for the two-dimensional two-electron coupled 
quantum dot and briefly describe the method used to cal- 
culate its electronic structure. In Section III, we calculate 
the spatial entanglement in the presence of one impurity 
and discuss its relation to the exchange coupling. Sec- 
tion IV, contains calculations of the oscillator strength, 
for a range of parameters of the system, that show the 
modifications of the optical properties in the presence of 
a charged impurity, with the aim of allowing to correlate 
optical measurements with the degree of entanglement 
of the system. Finally, In Section V we summarize the 
conclusions with a discussion of the most relevant points 
of our analysis. 



II. MODEL AND CALCULATION METHOD 

We consider two laterally coupled two-dimensional 
quantum dots whose centers are separated a distance d 
from each other, and containing two electrons. In quan- 
tum dots electrostatically produced, both their size and 
separation can be controlled by variable gate voltages 
through metallic electrodes deposited on the heterostruc- 
ture interface. The eventual existence of doping hydro- 
genic impurities, probably arising from Si dopant atoms 
in the GaAs quantum well, have been experimentally 
studied". Furthermore, some avoided crossing and lifted 
degeneracies in the spectra of single-electron transport 
experiments have been attributed to negatively charged 
Coulomb impurities located near to the From 
fitting the experimental transport spectra to a single- 
electron model of softened parabolic confinement with a 
Coulombic charge q, a set of parameters are obtained; 
among them, a radius of confinement of 15.5 nm, a con- 
finement frequency tku — 13.8 meV and an impurity 
charge of approximately 1 or 2 electron charges were 
found to be consistent. Actually, the uncertainty in 
the parameters and the suppositions introduced in the 
model does not allow one to precisely ensure the impurity 
charge, with the screening probably reducing its effective 
value to less than an electron charge. Therefore, we con- 
sider the charge of the doping atom Ze as a parameter 
varying in the range < Z < 1, in order to explore its 
effect on the properties of the system. 

In this work we model the Hamiltonian of the two- 
dimensional two-electron coupled quantum dot in pres- 
ence of charged impurities within the single conduction- 
band effective-mass approximation^, namely, 



h(r) 



2m 



-V 2 + V L (r) + V R (r) + V A (r), (2) 



where h(r) is the single-electron Hamiltonian that in- 
cludes the kinetic energy of the electrons, in terms of 
their effective mass m*, and the confining potential for 
the left and right quantum dots Vj, and Vr , and the inter- 
action of the electrons with the charged impurities, Va- 
The last term of the Hamiltonian, Eq. ([TJ , represents the 
Coulomb repulsive interaction between both electrons a 
distance ri2 = |r2 — ri| apart from each other, within a 
material of effective dielectric constant e. We model the 
confining potentials with Gaussian attractive potentials 



V l (r) = -V exp^- — \r-K l \ 2 j, (i = L,R), (3) 

where Rl and Rr are the positions of the center of the 
left and right dots, Vq denotes the depth of the potential 
and a can be taken as a measure of its range. Along this 
work, we will consider a single impurity atom centered 
at Ryi, and modeled as a hydrogenic two-dimensional 
Coulomb potential 



Va(t) 



Ze 2 
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(4) 



Since the Hamiltonian does not depend on the electron 
spin, its eigenstates can be factored out as a product of 
a spatial and a spin part 



%(r l ,T2,m ai ,m aa ) = *f (r 1 ,r 2 )xs 1 



M: 



(5) 



where 5 = 0,1 for singlet and triplet states, respectively, 
and M — m Sl + m s . 2 is the total spin projection. 

The eigenstates of the model Hamiltonian can be ob- 
tained by direct diagonalization in a finite basis set 2§ . 
The spatial part is obtained, in a full configuration inter- 
action (CI) calculation, as 



Nconf 

E c m„^n(ri,r 2 ) (6) 



where N con f is the number of singlet (S — 0) or triplet 
(S = 1) two-electron configurations $„(ri,r2) consid- 
ered, and n — is a configuration label obtained from 
the indices i and j from a single electron basis, i.e., 



H = h(n) + h(r 2 ) + 



where = (xi,yi) (i — 1,2) and 



47T££ r-12 ' 



(1) 



^(n,r 2 ) = -j= [^(n)^(r 2 ) + (1 - 25)^(n)^(r 2 )] 

(7) 

for i ^= j, and $^ _0 (ri,r 2 ) = </>i(ri)</>i(r 2 ) for the doubly 
occupied singlet states. 
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We chose a single-particle basis of Gaussian func- 
tions, centered at the dots and atom positions Rp (P = 
L,R,A), of the type28Jtt 

&(r) = Nx mi y ni exp (-a^r - R P | 2 ) , (8) 

where N is a normalization constant, and ti = rrii + 
rii is the z-projection of the angular momentum of the 
basis function. The exponents ai were optimized for a 
single Gaussian well and a single atom separately, and 
supplemented with extra functions when used together. 
For our calculations a basis set of 2s2p functions for the 
dots, and bsbpldlf for the atom was found to achieve 
converged results for the energy spectrum. 

The numerical results presented in this work refers to 
those corresponding to the parameters of GaAs: effec- 
tive mass m* — 0.067m e , effective dielectric constant 
e = 13.1, Bohr radius a* B — 10 nm and effective atomic 
unit of energy 1 Hartree* = 10.6 meV— >22. The depth 
of the Gaussian potentials modeling the dots are taken 
as Vb = 4 Hartree* = 40.24 meV, and its typical range 
a = \/2a* B = 14.1 nm. 



III. ENTANGLEMENT ENTROPY AND 
EXCHANGE COUPLING 

The proposed applications of QDs for quantum com- 
puting require a large exchange coupling between elec- 
trons along separated regions of space. To some extent, 
both requirements compete with each other. In a sim- 
ple picture, one could have a large exchange coupling for 
electrons doubly occupying the same dot or atom. In 
such a case, the singlet state has the form of a product 
wave function </?o( r i)¥>o( r 2) with the corresponding sin- 
glet spin function; the lowest triplet state, however, has 
the form of the antisymmetrized product of two single- 
electron functions, <po( r iVi( r 2) — <Po{ r 2) l fi{ri), of differ- 
ent single-particle energies Sq and e\. Thus, the triplet 
state will have a quite higher energy than the singlet 
state, thus giving a large exchange coupling. Neverthe- 
less, such a large coupling is not favorable for quantum 
computing tasks because the states are localized in space. 
Using electron states as qubits requires, for instance, the 
feasibility to detect the single or double occupancy of two 
quantum dots, separated a measurable distance, while 
keeping both electrons correlated. 

As the interdot separation increases, the electron- 
electron interaction diminishes and its relative impor- 
tance with respect to the confining potential tends to 
vanish. In the limit of large interdot separations, the 
Coulomb repulsion is minimized by singly occupying each 
quantum dots with an electron. In such a limit, the 
energies of both the singlet (+) and triplet (— ) states, 
yo( r i)vi( r 2) ± Vo( r 2)'/ 5 i( r i) approach the sum of singly 
occupied dots and their difference J tends to zero. In 
other words, the best conditions for applications to quan- 
tum information processing arises from a compromise be- 



tween a high spatial correlation of pairs of electrons at 
the longest possible lengths where the exchange coupling 
J is still sizeable. This behaviour is illustrated in Fig. [TJ 
assuming a positively charged impurity of one electron 
charge [Z = 1). 

Fig Q] shows the singlet and triplet ground-state ener- 
gies for the double QD, separated a distance d = 30 nm, 
as a function of impurity position xa- The inset shows 
the behaviour of the singlet-triplet exchange coupling as 
a function of the impurity position. These results are in 
qualitative agreement with those of Ref.— . The singlet- 
triplet exchange coupling is strongly affected for the pres- 
ence of the charged impurity, it has the maximum value 
when the impurity is centered in between the two dots, 
and it has a minimum close to zero when the impurity is 
located at xa — d. Expectedly, the splitting goes asymp- 
totically to the impurity-free double QD case when the 
impurity atom is located far away from the double QD 
system. We shall show below that the impurity positions 
that give high energy splitting, i.e., those near to the mid- 
dle of the interdot distance, correspond to a two-electron 
ground state whose spatial wave function is highly local- 
ized at the impurity atom, thus having a small spatial 
entanglement. 

In what follows we shall restrict ourselves to the impu- 
rity located along the interdot x-axis, Ra = (xa, 0). 
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FIG. 1. (Color online) Calculated ground-state energy of the 
two-electron double QD, for an interdot distance d = 30 nm, 
with one impurity atom of charge Z = 1 located along the 
interdot axis as a function of impurity position. The black 
line (red-dashed line) shows the singlet (triplet) ground-state 
energy. The inset shows the singlet-triplet exchange coupling 
J as a function of impurity position. 

We shall study now how the degree of spatial quan- 
tum correlation of two electrons in the coupled QD is 
modified by the position and charge of a screened atomic 
impurity. As mentioned above, the eigenstate wave func- 
tions can be factorised in its orbital and spin part. For 
the ground state, the spin part is a singlet wave function, 
which is maximally entangled and constant. Therefore, 
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throughout this work, we will only consider the spatial 
entanglemen t 10 ' 13 ^"—. 

The Von Neumann entropy of the reduced density ma- 
trix quantifies the entanglement for a bipartite pure state 
and can be calculated usin g 13 ' 15 " — 

S=_Tr(p red log 2 p red ), (9) 

where p rod = Tr2 1^) ("l/l is the reduced density operator, 
*S> is the two-electron wave function and the trace is taken 
over one electron. The Von Neumann entropy could be 
written as 

S = -^Vlog 2 A 4 , (10) 

i 

where A^ are the eigenvalues of the spatial part of the 
reduced density operator 

J ^(n, ri)^(ri) dr[ = A^n) . (11) 

where 



/ cd (r 1 ,r' 1 )= / ^(n.ra^W.raJdra. (12) 

Fig. [5] shows the Von Neumann entropy for two electrons 
in the double QD as a function of its interdot distance, in 
absence of impurity (black-dashed line) and with atomic 
charges Z = 1 (blue-dashed doted line) and Z = 0.1 
(red solid line) located at the center of the double QD. 
In all the cases it is observed that, for small interdot 
separations, the entropy is small, smoothly increasing 
with the interdot separation. The increase of the spa- 
tial entanglement is due to a gradual delocalization of 
the ground state wave function. For interdot distances 
between 20 and 40 nm, there is a large increase of the 
entanglement entropy, signaling a qualitative change of 
the ground state wave function from a atomic doubly oc- 
cupied state to a state with both dots singly occupied, 
reached at large interdot separations (d > 50 nm) , where 
the entropy saturates to its maximum 5 = 1. The varia- 
tion of S is similar for all the cases, although the presence 
of the impurity decreases the entanglement for every in- 
terdot separation, due to the fact that the atomic poten- 
tial contributes to localize the electronic density at the 
center of the system. 

The effect of the charge and location of the impurity on 
the spatial entanglement, for fixed QDs geometry, can 
be observed at Fig. [3] where the entropy is depicted 
as a function of the impurity position. The separation 
between the two QDs is kept fixed at 30 nm, and two 
limiting cases are considered: a highly screened atomic 
charge Z = 0.1 and a unscreened charge Z — 1. In both 
cases, the entanglement entropy increases as the impu- 
rity moves off the center of the double QD until a position 




d [nm] 



FIG. 2. (color on-line) Von Neumann entropy of the reduced 
density matrix for the two-electron coupled QD as a function 
of the interdot distance. Black-dashed line shows the entropy 
when there is no impurity present in the sample, the red line 
corresponds to the entropy when a single impurity Z = 0.1 is 
located at x = y = and blue-dotted dashed line corresponds 
to Z = 1.0. 
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FIG. 3. Von Neumann entropy of the reduced density matrix 
for the two-electron coupled QD as a function of the impurity 
position (along the interdot axis) for d — 30 nm and different 
impurities strength: (a) Z — 0.1 and (b) Z = 1.0. 



where S reaches a maximum, finally decreasing to a value 
S = 0.53, when the atom is distant from the dots (x > 40 
nm). 

The minimum and maximum of entanglement produced 
by the small charge Z = 0.1 are less pronounced than 
those due to the highly charged impurity Z = 1. This 
modulation of the entanglement by the impurity position 
reflects the existence of two regimes: one of low entan- 
glement for impurities at (or near to) the center of the 
interdot distance, and another of higher (but not max- 
imun) entanglement for atomic positions external to the 
interdot segment. In Fig. [3J this two regions are the 
ones to the left and the right of the bell-shaped peak of 
S, repectively. The peak position itself depends on the 
magnitude of the charge. For small charges, the maxi- 
mum degree of entanglement occurs at xa ~ 17 nm, that 
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is, close to the center of the dot to the right. For the large 
charge Z = 1, however, the peak of S occurs at x s» 30 
nm. The rationale for it is that, for low charged impuri- 
ties, the atomic potential is a weak perturbation to the 
QDs potential wells. Therefore, the entropy varies in a 
small range (0.47 < S < 0.56) around the impurity-free 
case S = 0.53. For large impurity charges, nevertheless, 
the atomic potential is strong and the position of its cen- 
ter greatly determines the spatial wave function. The 
range of atomic positions (0 < xa < 20 nm) along which 
S remains low, can be understood as due to the localiza- 
tion of the electrons close to the atom. When the atom 
is inside one of the dots (QD#), the atomic potential 
reinforces that of the dot well and the electron density 
localization, thus giving a low degree of entanglement. 
When the atom moves towards outside the double QD, 
the strength of the double well competes with the large 
atomic potential until an atom-double QD distance of ca. 
30 nm, where becomes energetically convenient to delo- 
calize the electron wave function, resembling the doble 
QD bond in absence of impurity. 

To show clearly the influence of the atomic charge 
on the wave function, let us consider two a bit less ex- 
treme situations: Z = 0.2 and Z = 0.8. Fig. [4] shows 
the ground state electron density along the interdot axis 
when the impurity atom is located at xa — 15 nm, for 
three different interdot separations, d — 15, 25 and 40 
nm. The panels to the left show that for the small charge, 
as the QDs separate from one another, the electron den- 
sity develops peaks located at the potential well centers. 
For the large charge Z — 0.8, however, the density is al- 
ways peaked at the impurity position. Therefore, in this 
last case, the presence of the impurity could spoil the 
performance of the device for quantum computing tasks 
due to the high localization of the electron density entails 
a low degree of entanglement. 

Fig. E] shows the dependence of the Von Neumann en- 
tropy on the impurity charge and interdot distance for 
a given position of the impurity atom: xa = 15 nm. 
Ideally, providing that the value of the impurity charge 
could be measured in a given sample, one would be able 
to choose the optimal interdot distance for a given degree 
of entanglement. The figure clearly shows the aforemen- 
tioned regimes of weak (Z < 0.6) and strong [Z > 0.6) 
impurity potential, corresponding to low and high degree 
of entanglement, respectively. For a given (fixed) small 
impurity charge Z, the entropy increases monotonically 
as the interdot distance d increases. On the other hand, 
for a given large Z, by increasing the distance d, the en- 
tropy increases for small distances d up to a maximum 
value, diminishes to a minimum and sharply increases 
again until its asymptotic impurity- free value 5 = 1. 

Fig. [6] shows the dependence of the entropy on the in- 
terdot distance, for different impurity positions {xa = 15 
and 20 nm) and charges (Z = 0.1, 0.5, 0.7 and 1). The 
corresponding variation in absence of impurity is also rep- 
resented in dashed lines for reference. It can be seen two 
qualitatively distinct behaviours for small [Z = 0.1, 0.5) 
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FIG. 4. (color on-line) Ground state one-electron density 
along the interdot axis. Left panels show the weak impu- 
rity limit Z — 0.2 and right panels the strong impurity limit 
Z = 0.8: (a) and (d) d = 15nm, (b) and (e) d = 25nm and 
(c) and (f ) d = 40nm . Red circles show the impurity position 
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FIG. 5. (color on-line) Contour map for the Von Neumann 
entropy of the reduced density matrix for the two-electron 
coupled QD as a function of impurity strength and interdot 
distance for xa = 15 nm. 



and large (Z — 0.7, 1) charges . The monotonic increase 
of S with d is characteristic of the weak atomic poten- 
tial; separating the QDs with a small atomic charge in 
between of them, produces little changes in the electron 
distribution as compared with the impurity-free double 
QD. On the other hand, strong atomic potentials induce 
a modulation of the entropy as d increases; for small val- 
ues of d, all three potentials are close to each other and 
the electron density localizes around their centers. For 
large interdot distances, the energy of the system is min- 
imized by decreasing the electronic repulsion, i.e., by de- 
localizing the wave function and, hence, increasing its 
entanglement. 

Fig. [7] shows this effect on the exchange coupling corre- 
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FIG. 6. (color on-line) Von Neumann entropy of the reduced 
density matrix for the two-electron coupled QD as a func- 
tion of the interdot distance for xa = 15 nm ((a) and (c)), 
x = 20 nm ((b) and (d)) and different values of the impurity 
Strength, (a) and (b) show the weak impurity limit (Z — 0.1 
and Z = 0.5) while in (c) and (d) we observe the strong im- 
purity limit (Z — 0.7 and Z = 1.0). The black-dashed line 
shows the behavior of the entropy when there is no impurities 
in the sample. 



ponding to situations of Fig. UJa) and[5]Jc), having the 
atom at x a = 15 nm. It can be observed that S and 
J have, roughly, opposite variations; whence the atomic 
potential is weak, S increases and J decreases as the 
QDs separate from each other. In the regime when the 
atomic potential is strong, the maximum of S occurs at 
the minimum of J and reciprocally; furthermore, at large 
QDs separations, as the entropy goes to its asymptotical 
value S = 1 the exchange coupling tends to zero. Then, 
for specific quantum information applications, it could 
be desirable to tune the interdot distance for harnessing 
one or both properties. 

The variety of behaviours of the degree of spatial en- 
tanglement with the various parameters of the system, 
described in this section, is rooted in the spatial distri- 
bution of the electron wave function. We shall discuss 
in the next section a relation with an optical property, 
like the oscillator strength, in order to provide a feasible 
connection with measurable magnitudes. 



IV. IMPURITY EFFECT ON THE OPTICAL 
PROPERTIES 

The optical susceptibility of a system depends on its 
transition amplitude for the interaction of its dipole mo- 
ment with the optical electric field between two singlet 
states iffi and ^>j, say the ground and an excited states, 
and the corresponding energy differences. The oscillator 
strength for an electric field applied along the interdot 
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FIG. 7. (color on-line) J as a function of the interdot distance 
d with an impurity center located in xa = 15nm. On the left 
panel (a) we observe the singlet-triplet coupling for Z = 0.1 
(red circles) and Z — 0.5 (blue squares). In (b) we show 
the singlet-triplet coupling for Z = 0.7 (green up triangles) 
and Z = 1.0 (orange down triangles). The black-dashed line 
represents the singlet-triplet with no impurity. 
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(13) 



takes both magnitudes into account and provides infor- 
mation on the feasibility of such optical excitations. 

We study here how the impurity affects the oscillator 
strength of the double QD. The dots are kept 30 nm sep- 
arated from each other and the position of the impurity 
xa is varied from the center of interdot segment (xa = 0) 
to a large separation from the dots (xa = 70 nm), includ- 
ing the case of the impurity centered in one dot (xa = 15 
nm). The charge Z of the atom is varied from Z = 0.1 
(highly screened impurity) to Z = 1 (low screening). For 
the system considered, the oscillator strength between 
the ground and the first excited singlet states, /oi, is 
the dominant contribution with respect to all others fij. 
The precision of the calculation was checked by verifying 
the Thomas-Reiche-Kuhn sum rule, J^ij fij — with 
N = 2 being the number of electrons in the system. The 
results are shown in Fig. [8] toghether with the entangle- 
ment entropy for the same atomic positions and charges. 

The cases of weak and strong electron-atom Coulomb 
interaction are clearly distinguishable. In the regime of 
small impurity charge (Z < 0.6), the oscillator strength 
/or varies approximately with a quadratic dependence 
on xa', i.e., it starts from /ni ~ 2, reaches a minimum 
around x a — 15 nm, to finally increases up to a value of 
2, at nearly xa = 30 nm. The larger the impurity charge 
Z, the more pronounced the minimum of /oi . Placing the 
atom further away from the double dot system (xa > 30 
nm) does not change /oi- 

On the other hand, in the regime of high impurity 
charge (Z > 0.6), the oscillator strength /oi exhibits 
richer features as compared to the small charge case. The 
most remarkable behavior correspond to Z = 1 which 
successively shows a similar decreasing, from /oi = 1.4 
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FIG. 8. Von Neumann entropy of the reduced density ma- 
trix (upper panel) and Oscillator strength /oi (lower panel) 
between the ground and first excited singlet state of a double 
quantum dot with an impurity of charge Z as a function of 
the impurity position xa, for different values of Z. 



at xa = 0, to foi — 0.6 at xa — 15 nm, followed by an 
increase up to xa = 25 nm, a small plateau around 30 
nm, a peak at xa = 33 nm, a minimum of /oi ~ 1 at 
40 nm, finally approaching the saturation value /oi = 2 
for xa > 50 nm. For intermediate 0.5 < Z < 1 values, 
a gradual transition between both regimes is observed; 
namely, by decreasing Z from 1 to 0.5, the minimum of 
the region xa ~ 40 nm becomes shallower, the peak is 
softened, and the plateau merges with the minimum oc- 
curing at 15 nm, thus giving the flat minimum of the 
weak impurity regime. 

It should be noted that the range < x < 15 nm 
corresponds to an impurity atom located in between the 
dots, while for xa > 15 nm, the atom is outside the 
segment defined by the interdot centers. Consequently, 
the existence of an impurity into the system would 
cause a diminishing in /oi and, therefore, in the light 
absorption or emission of the double dot device. This 
effect is stronger the closer is the atom to one dot. The 
most favourable situation for optical excitation (high 
fax) correspond to an impurity centered in between the 
dots or outside the interdot separation, faraway from 
any of them. 

We shall discuss in the following, the behavior of the 
oscillator strength as due to changes in the electronic 
structure induced by the variation of the position of the 
impurity, starting by the most striking case of a highly 
charged impurity Z = 1. We displace the atom along the 
line joining both dots, which we take as the x-axis; there- 
fore, we consider the two-particle wave function along the 



x axis for the coordinates x\ and xi of each electron 

*i(xi, x 2 ) = *i(n, r 2 ) = *<(ii,0; x 2 , 0), (14) 

for the two lowest singlet states \&j(ri, r 2 ), i — (ground 
state) and i = 1 (first excited state). The function 
^(xi, X2), represented as a two-dimensional plot in the 
(xi, X2)-plane, allows one to visualize the most relevant 
configurations contributing to the total wave function. 
Because of the permutation symmetry, the spatial wave 
function satisfies i ff(xi,x 2 ) = ^(x2,%i)i thus becoming 
symmetric under reflection with respect to the diago- 
nal X\ = X2- Large values of ^(x, x), along this diag- 
onal, correspond to ionic or doubly occupied configura- 
tions. In contrast, large density values 9(x,—x) along 
the X\ = —x 2 diagonal, corresponds to configurations 
where the electrons are mostly in opposite (left and right) 
half-planes. 

In the present calculations, the x coordinates of the 
centers of the left and right dots xl = —15 nm, xr = 15 
nm are held fixed while that of the atom, xa = x, varies. 
Large values of 'F (xl , xl ) , ^ (xr , xr) or 'F (xa , xa ) entail 
a doubly occupied configuration at the left dot, the right 
dot or the atom, respectively. 

On the other hand, a configuration of one electron in 
the atom and the other in a bond (antibond) between 
the left and right dots, would be represented by 

*(xi, x 2 ) = [cl<Pl{xi) ± c R (p R {x\)] (Pa(x2) + (xi «-» x 2 ), 

(15) 

where the last term represents a term similar to the 
first one with the variables interchanged, and ip a is a 
wave function centered around x a (a = L,R,A). Then, 
i £(x\,X2) will have large values close to (xl,xa) and 
{xr, xa) with the same or opposite sign for a bond or 
antibond, respectively. 

Figs. IH1 and [TU] show the plot of the ground state 
\l/rj(iri, X2) and the first singlet excited state \{x\,X2) 
from our calculations. 

For a single-electron symmetric double dot system 
without impurity, the ground and first excited states are 
the bonding and antibonding states formed from the lin- 
ear combination of orbitals centered at each dot. For the 
two-electron symmetric double dot system, Figs. [S] and 
[TU1 show that when the atom is at the center of the line 
joining both dots (x = 0), the two-particle wave function 
of the ground state (excited state) roughly corresponds to 
one electron in the atom and the other in the bond (an- 
tibond) of the double dot system, Eq. (|T5|) . Therefore, 
the matrix element ("Fol^i + a^l^i) roughly correspond 
to the sum of those for the atom and the double dot 
separately. 

At x — 15 nm, the atom is at the center of the dot 
to the right, the system becomes very asymmetric, with 
the potential of dot to the right deeper than the one to 
the left due to the contribution of the attractive impu- 
rity. The bond becomes a doubly occupied state localized 




FIG. 9. (Color online) Countour plot in the (0:1,0:2)- 
plane of the two-electron ground state singlet wave function, 
^0(2:1, X2), along the interdot axis x for the ground state of 
the doped double quantum dot. Gaussian wells are centered 
at Rl = (xl,Q) and Rr = (xr,0) and the impurity atom of 
charge Z = 1, at Ra = (xa,0). The vertical and horizontal 
dashed lines Xi — xl,xr or xa (i = 1, 2), signals the condi- 
tion where one electron (electron 1 or 2, respectively), is at 
the center of the dot to the left, to the right or at the impurity 
atom. The centers of the dots xr = — xl = 15 nm are held 
symmetrical with respect to the origin of coordinates. The 
atom is successively placed at xa = 0, 15, 20, 25, 30, 33, 40 
and 50 nm. 



close to the center of the combined potential (QDr and 
impurity), while the antibond becomes more localized 
around to QD^, due the orthogonality condition, what 
lowers /01. The behavior in the range < x < 15 nm 
reflects this gradual change. 

From 15 to 30 nm, the effect of the impurity turns 
weaker as the atom moves away, and the two dots be- 
comes more symmetric again; this redistribute the charge 
towards QD^, recovering some bonding and antibonding 
character for \I/o and ^i, respectively. Such a config- 
uration favors an increase of /oi- Furthermore, as one 
electron remains in the atom, which is farther from the 
origin, the matrix element of x becomes larger than the 
one corresponding to the atom at the origin. The oscil- 
lator strength has a peak at 33 nm roughly increasing 
quadratically with the position of the atom as a conse- 
quence of the stretching of the charge. 

After 33 nm, the electron in the atom cannot be re- 
tained by the impurity potential, thus "Fo approaches to 
a configuration with one electron in each dot. Neverthe- 
less, the excited state >Fi still have a configuration where 
the atom is occupied, what lowers /oi- 

For atom positions further than 40 nm, the excited 
state also releases its electron and the double quan- 
tum dot becomes even more symmetric, approaching 
its behavior in the absence of impurity, thus approach- 



FIG. 10. (Color online) Same as Fig. [5] for the first excited 
two-electron singlet wave function $5(a;i, x%). 



ing its value /oi = 2. The limit of isolated dots is 
clearly seen in Figs. H] and [TU1 where for x > 50 
nm, the ground and first excited states are, approxi- 
mately, <F ~ [^l(xi)(Pr(x 2 ) + ifR(x 1 )ip L (x 2 )]/V2 and 

*1 ~ [fL{xi)ifL{x2) ~~ ipR(xi)ip R (x 2 )]/V2. 

As seen from Fig. [8j the oscillator strength for impu- 
rity charges smaller than Z = 1, has simpler features. 
Basically, they start from a value /oi slighly less than 2, 
decreases until a minimum as the atom approaches one 
dot, say QD#, and increases again smoothly until reach- 
ing the asymptotic impurity-free value of 2. 

The oscillator strength is a highly sensitive property 
to the presence of the impurity. A value of foi close 
to 2, occurs either when the impurity is weak wherever 
it is located, or when a highly charged impurity atom 
is far away from the double QD. Both cases are situa- 
tions where the impurity is a perturbation for the coupled 
QDs and, therefore amenable of use in quantum comput- 
ing. On the contrary, deviations of the oscillator strength 
from a value of 2, provides an indication of a breakdown 
of the possiblity to consider the system as a double QD. 



V. CONCLUSIONS 

In this work, we have studied the influence of a 
Coulombic atomic impurity on the entanglement entropy 
of two-dimensional two-electron double quantum dots. 
The electronic structure was calculated by using a con- 
figuration interaction method with a Gaussian basis set 
expansion. The degree of entanglement becomes highly 
modulated by both the location and charge screening of 
the impurity atom. Two regimes are clearly identified: 
one of low entanglement and other of high entanglement, 
with both of them mainly determined by the magnitude 
of the charge. The exchange coupling between the elec- 
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trons, being proportional to singlet-triplet exchange cou- 
pling , has an opposite behaviour with respect to the 
one of the entropy. The efficient use of double QDs with 
impurities, in specific quantum information processing 
tasks could require the tuning of the interdot separation 
or the quantum well depths, for optimizing the harness- 
ing of the entanglement, the exchange coupling or both. 
Finally, the magnitude of the oscillator strength of the 
system could provide an indication of the presence and 
characteristics of impurities that could largely influence 
the degree of entanglement of the system. It is clear 
that experimentally obtained optical properties can help 



in the design of double QDs with desirable properties in 
order to use them in quantum information tasks. 
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